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By using the cubic spline interpolation method, we reconstruct the shape of the primordial scalar 
and tensor power spectra from the recently released Planck temperature and BICEP2 polarization 
cosmic microwave background data. We find that the vanishing scalar index running (dn;/d Ink) 
model is strongly disfavored at more than 30 confidence level on the k = 0.0002 Mpc™! scale. 
Furthermore, the power-law parameterization gives a blue-tilt tensor spectrum, no matter using only 
the first 5 bandpowers nz = 1.207924 (95%CL) or the full 9 bandpowers nz = 1.247923 (95%CL) of 
BICEP2 data sets. Unlike the large tensor-to-scalar ratio value (r ~ 0.20) under the scale-invariant 
tensor spectrum assumption, our interpolation approach gives ro.oo2 < 0.060 (95%CL) by using the 
first 5 bandpowers of BICEP2 data. After comparing the results with/without BICEP2 data, we 
find that Planck temperature with small tensor amplitude signals and BICEP2 polarization data 
with large tensor amplitude signals dominate the tensor spectrum reconstruction on the large and 
small scales, respectively. Hence, the resulting blue tensor tilt actually reflects the tension between 


Planck and BICEP2 data. 


PACS numbers: 


I. INTRODUCTION 


Recently BICEP2 experiment reported an excess 
of Cosmic Microwave Background (CMB) B-mode po- 
larization power spectrum over the base lensed-ACDM 
expectation in the range 30 < £ < 150, inconsistent with 
the null hypothesis at a significance of > 50. Since the 
single field slow-roll inflationary model predicts a peak 
around multipole £ ~ 80 in the B-mode auto-correlation 
(BB) spectrum seeded by the primordial gravitational 
wave/tensor perturbation mode, the BICEP2 results are 
believed as the first indirect detection of the primordial 
gravitational wave. Under the assumption of power-law 
scalar and scale-invariant tensor spectra, the observed 
B-mode power spectrum is well-described by a lensed- 
ACDM++tensor theoretical model with tensor-to-scalar 
ratio r = 0.20790! with r = 0 disfavored at 7.00 con- 
fidence level. 

However, the scientific results of BICEP2 data are in 
tension with those from other CMB experiments, such 
as Planck d. The first discrepancy is in the amplitude 
of the scale-invariant tensor spectrum, which is described 
by the tensor-to-scalar ratio r = A,/A,. Unlike the scalar 
perturbatoins, due to the absence of acoustic oscillation 
mechanism the tensor contributions to the temperature 
CMB spectrum are rapidly washed out inside the hori- 
zon at electron-proton recombination epoch (€ > 200) 

Hence, the temperature anisotropies on the large 
scales are the mixture of scalar and tensor contributions. 
Furthermore, if one assumes the simple power-law form 
of the primordial scalar power spectrum, i.e., no scalar 
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index running dns/dlnk = 0, the precisely measured 
higher multipoles by Planck put stringent constraints on 
the scalar amplitude A, and index n,. Therefore, in or- 
der to explain the observed power deficit in the low-@ 
regime by Planck, one has to suppress the tensor spec- 
trum amplitude. Consequently, from only the tempera- 
ture anisotropy measured by Planck, one has r < 0.11 
at 95%, which is in a “very significant” tension (around 
0.1% unlikely) with BICEP2 results g. As stressed by 
the BICEP2 team, however, this tension could be rec- 
onciled by adding the running of the scalar index which 
is the degree of freedom to suppress the scalar tempera- 
ture anisotropy in low-/ regime. Then, combining Planck 
with WMAP low-@ polarization [5] and ACT 6] /SPT [+ 


[9j high-@ data, one could get r < 0.26 at 95%. Besides 


that, several other possible solutions to this tension have 
been proposed, such as step feature spectra LHA, fast- 
slow roll model {13}, anti-correlation scalar iso-curvature 
initial condition 14 05), sterile neutrino species Hig, 
sudden change in speed_of inflaton or Lorentz violation 
[19}. And see also for other possibilities. 


The second tension is between the observed blue-tilt 
tensor spectrum (ng > 0) [RH30] by BICEP2 and the 
red-tilt one (nz = —7r/8) predicted by the standard infla- 
tionary paradigm. Generally, the blue tensor spectrum 
asks for violation of Null Energy Condition (NEC), which 
is equivalent to p+ P < 0 (or H > 0) in the flat universe. 
There exist several NEC violation inflationary models in 
the literature, such as super- inflation [31] Bi, phantom infla- 
tion Bd, G-inflation etc. Some alternative] oe 
of inflation, ln a “string gas ate ee 
ing universe or other possibilities U2 (see the 
references Heel might be helpful to a ‘te tension. 


In this paper, we start from a purely phenomenological 
point of view to reconstruct the shape of primordial scalar 
and tensor spectra from Planck temperature and BICEP2 
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polarization data. 


II. PARAMETERIZATION OF PRIMORDIAL 
SPECTRA AND DATASETS 


In order to reconstruct a smooth spectrum with the 
continuous first and second derivatives, in this paper 
we adopt the cubic spline interpolation method, which 
has been used to analyze WMAP or Planck tempera- 
ture and polarization data [44 {53}. Beside the cubic 
spline interpolation reconstruction, there also exist other 
model independent algorithms, such as Bayesian evidence 
selected linear interpolation, (52), etc. Due to the 
fact that the amplitude of CMB anisotropy is so tiny, 
6T/T ~ O(10-*), and CMB observational windows cover 


| 


dln P(k) 


several orders of magnitude in spatial scale, it is reason- 
able to parameterize the logarithms of primordial spec- 
tra, which seed the CMB anisotropy, in the logarithms 
of fluctuation wavenumber, Ink. The method of cubic 
spline interpolation can be summarized as follows. 


First, we uniformly sample Npin points in the logarith- 
mic scale of wavenumber. Second, inside of the sampled 
bins Ink; < Ink < Inkj+1, we use the cubic spline inter- 
polation to determine logarithmic values of the primor- 
dial power spectrum. Third, the boundary conditions 
are adopted, where the second derivative is set to zero. 
For k < kı or k > ky,,, we fix the slope of the pri- 
mordial power spectrum at the boundaries and linearly 
extrapolate to the outside regimes. Mathematically, the 
corresponding formula could be written as 


dink |, In = + In P(k1), k < kı; 
$ 
In P(ki), k € {ki}; 
mP(k) = 4 eubic spline, kid ké kua (1) 
din P(k) k 
— kny, In ENyin + InP(kyyin)s k> ENpin: 


This reconstruction method has three advantages: first 
of all, it is easy to detect deviations from a scale-invariant 
or a power-law spectrum because both the scale-invariant 
and power-law spectra are just straight lines in the In k- 
ln P plane. Second, negative values of the spectrum can 
be avoided by using In P(k) instead of P(k) for splines 
with steep slopes. Finally, the shape of the power spec- 
trum reduces to the scale-invariant or power-law spec- 
trum as a special case when Npin = 1, 2, respectively. 


Since the purpose of this work is to reconstruct the 
scalar and tensor spectra, we need to adopt different 
sampling logarithms based on the different observational 
windows. For the primordial scalar curvature spec- 
trum, its constraints are mainly driven by the CMB 
temperature modes. With Planck sensitivities, we uni- 
formly sample 3 bins ranged in Ink € (—8.517, —1.609), 
which corresponds to k € (0.0002,0.2) Mpce~!. For the 
tensor spectrum, we adopt two uniformly logarithmic 
sampling strategies, one is corresponding to the scales 
k € (0.002,0.03) Mpc~', the other is k € (0.002, 0.02) 
Mpc~!. This is because the BICEP2 B-mode polariza- 
tion data, which is a very sensitive probe for the primor- 
dial tensor spectrum, have an excess of B-mode power in 
all the range of polar angle mulitpoles (20 < @ < 340). 
And moreover, compared with the first 5 bandpowers, 
which is in the range of (20 < Z < 200), the power in the 
second 4 bandpowers has extraordinary excess over the 
base lensed-ACDM expectation. This extraordinary ex- 
cess might arise from exotic physical origin beyond stan- 
dard inflationary paradigm or from some unresolved fore- 


m 


ground contaminations. Given this consideration, in this 
work we take two different choices of BICEP2 data, the 
first is using the full 9 bandpowers, and the second is to 
use the selected first 5 bandpowers. Hence, we have to 
adjust our sampling logarithms as mentioned above. 


In the left part of this section, we would like to briefly 
review the data sets we used. First of all, we utilize 
the Planck TT power spectra, namely, for low- modes 
(2 < £ < 50) via all the 9 frequency channels ranged 
from 30 ~ 353 GHz, for high- modes (50 < £ < 2500) 
through 100, 143, and 217 GHz frequency channels ! [55]. 
Second, in order to break the well-known parameter de- 
generacy between the re-ionization optical depth and the 
amplitude of CMB temperature anisotropy, we also in- 
clude WMAP9 low-¢ temperature/polarization spectra 
(2 < £ < 32) [B]. In addition, we use BICEP2 polar- 
ization (EE,EB,BB) spectra from 9 (or 5) bandpowers of 
multipoles in (20 < ¢ < 340 or 200) of 150 GHz chan- 
nels 2. For the data analysis numerical package, we 
compute the CMB angular power spectra by using the 
public Einstein-Boltzmann solver CAMB and explore 
the cosmological parameter space with a Markov Chain 
Monte Carlo sampler, namely CosmoMC B). 


1 [http://pla.esac.esa.int/pla/aio/planckProducts.html 
2 http: //www.cfa.harvard.edu/CMB/bicep2/papers.htm1 


201709.00159v1 


chinaXiv 


Parameter | Range (min, max) 
In(10!° A?) 2.7, 4.0) 

Ns 0.9, 1.1) 
dns/dink (—1.0, 1.0) 

10.05 0.0, 2.0) 

ne (—1.0, 5.0) 
In(10"° A1) 1.0, 5.0) 
In(10°° A2) 1.0, 5.0) 
In(10*° A3) 1.0, 5.0) 
In(10*° B1)* |(—3.0, 3.0) /(—3.0, 3.0) 
In(10!° B2)* | (—2.0, 5.0) /(—3.0, 3.0) 
lIn(10*° B3)* |(—1.0, 5.0) /(—3.0, 3.0) 


TABLE I: List of the primordial spectrum parameters used in 
the Monte Carlo sampling. *The left parameter ranges are 
for the chains from Planck+WP+BICEP2 data compilation, 
and the right one are for those without BICEP2 data. 


II. RESULTS AND DISCUSSIONS 


In this section we will start with the scalar spectrum re- 
construction, and then turn to the tensor spectrum case. 
The prior ranges of the primordial spectrum parameters 
we studied are listed in Tab[I] Here we emphasize that 
the differences in the prior of tensor spectra amplitude 
at our cubic spline sampling knots when BICEP2 data 
are included, lie in the tension between Planck and BI- 
CEP2 data (we will show later). When our MCMC sam- 
pler investigates the wide parameter space spanned by 
(In B1, In Bg, In B3), at some points the resulting spectra 
are inconsistent with Planck TE cross-correlation data, 
in order to avoid this problem we have to adjust the prior 
ranges. But still the width of the priors are large enough 
and also the tensor amplitudes (In B),In Bz,In B3) get 
well constrained in these prior ranges as shown in Fig[9] 
So, we conclude that our prior choices will not affect the 
results significantly. 


A. Scalar spectrum reconstruction 


Since the main capability of scalar spectrum recon- 
struction is driven by the CMB temperature data, we 
firstly study the case without BICEP2 polarization, i.e. 
only with Planck temperature and WMAP9 low-é po- 
larization (WP) data sets. As mentioned in the previ- 
ous section, here we uniformly sample 3 points in the 
logarithmic scale of wavenumber, which are located at 
kı = 0.0002, kə = 0.0063 and k = 0.2 Mpc! with the 
logarithmic amplitudes In A;, In Ag and In As, respec- 
tively. And then, we sample the parameter space spanned 
by vanilla ACDM parameters without the scalar ampli- 
tude In A, and its index tilt ns and replacing them with 
In A,, In Ag and In As, hereafter we call this parameter 
compilation as ACDM-In A,-n,+ln Aj+In A2+ln A3. 

The marginalized mean scalar spectrum reconstructed 


T T 
—— ACDM - InA, — n, + nA, + InAy + InAg 


=) =) ACDM + dn,/dink 
35+ - - -ACDM 


10° 


k{Mpe'] 


FIG. 1: Reconstruction of primordial scalar spectrum without 
BICEP2 data. 


from Planck+WP data is represented by the black solid 
curve in Fig. [I] and the corresponding 1 ~ 30 er- 
ror bars at the sampling points are denoted by the 
blue, red and green segments, respectively. For com- 
parison, we also show the primordial scalar spectrum 
from the Planck marginalized mean vanilla ACDM and 
ACDM-+dn,/dink (scalar index running) with blue 
dashed and red dotted-dashed curves. From Fig. [I] 
we can see that, first, our cubic spline interpolation 
result mimics the ACDM+dn,/dlnk case; second, on 
k = 0.0002 Mpc~! scale, the simplest vanilla model is 
disfavored at nearly 20 level. 


r r 
—— ACDM - InA, — n, + InA + logy + InAs +r 
.——:ACDM + dn,/dink +r 
—--ACDM 
© k= 0.050, In [10"0P,(k)] = 3.089 
3.47 o k= 0.002,1n[100P,(k)]= 3.217 
lee © =k =0.050,In[10!°P,(k)] = 3.121 
:)] = 3.111 


( 
Pee o k=0.002,In{10!°P,(k 


In[10!P,(k)] 


2.67 


2.4 1 1 1 
10° 
k{Mpe'] 


FIG. 2: Reconstruction of primordial scalar spectrum with 
BICEP2 data. 


Adding the BICEP2 polarization data, we show the 
results in Fig. [2]for models including the tensor-to-scalar 
ratio r, i.e. ACDM-In A,-n,t+r+In Aı+ln A2+ln A3. We 
find that, first of all, due to the anti-correlation between 
scalar and tensor amplitudes (see the bottom left sub- 
panel of Fig. Blin Appendix [A}, the large value of tensor- 
to-scalar ratio discovered by BICEP2 data will lead to 
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the suppression of scalar amplitude on the large scales. 
This is also explicitly demonstrated in the top sub-panel 
of Fig. [8] in Appendix As a result of deficit of scalar 
power on the large scales, as shown in Fig. 2] the vanilla 
ACDM model (blue curve) is strongly disfavored with 
> 30 confidence level on the k = 0.0002 Mpc~? scale. 
A similar result is obtained by authors of . By using 
Planck and BICEP2 data, they found a distinct prefer- 
ence for a suppression of power in the scalar spectrum 
at large scales, k < 1078 Mpc™! via a linear spline re- 
construction method. Second, by assuming the scale in- 
variant tensor spectrum, our scalar spectrum cubic spline 
interpolation parameterization still gives a large tensor- 
to-scalar ratio r = 0.211649 at 95%C.L., see Tab. Tin 
Appendix [A] 


B. Tensor spectrum reconstruction 


In the previous subsection, we have assumed the ten- 
sor spectrum is scale-invariant. In this subsection we 
relax the assumption and use the same cubic spline 
interpolation method to reconstruct the shape of ten- 
sor spectrum. Unlike the CMB temperature spectrum 
which is mainly sourced by the primordial scalar pertur- 
bations, the B-mode polarization anisotropy seeded by 
tensor perturbation is only detected by BICEP2 on the 
large scales with the polar spherical harmonic mulitpoles 
ranged 20 < £ < 340 (9 bandpowers). As the case for the 
scalar spectrum, here the tensor spectrum is also uni- 
formly sampled with 3 points in the logarithmic scale 
of wavenumber, but only on the scales covered by BI- 
CEP2 observations. Consequently, we sample them at 
kı = 0.002, kə = 0.0077 and k3 = 0.03 Mpc7, respec- 
tively. The cosmological parameters we estimated are 
the 6 vanilla ACDM model parameters plus the extra 3 
tensor amplitudes In B,, In B2 and In B3. 


r T J 
— ACDM + InB, + InB + InB, A 


ap =- ACDM +m +r z 
O k=0.002,In{10!°P,(k)] = -0.4135 3 


In{10"" P;(k)] 


k{Mpc!} 


FIG. 3: Reconstruction of primordial tensor spectrum with 
BICEP2 data (9 bandpowers). 


In Fig. [3] we plot the primordial tensor spectra from 
the best-fit model of ACDM+In Bı+ln B2+ln B3 (black 


solid curve) and ACDM+r-+n; (red dotted-dashed curve) 
as well as the error bars at the sampling points of the cu- 
bic spline interpolation method. First of all, both the 
standard power-law and our cubic spline parameteriza- 
tions favor a blue-tilt tensor spectrum. The former (see 


Tabl] in AppendixB) reports 


ro.oo2 < 0.061, 95%CL , (9 bandpowers) , (2) 
me = 1.244528 ,95%CL , (9 bandpowers) , (3) 


and our cubic spline interpolation method gives 
r0.002 < 0.064 ,95%CL , (9 bandpowers) . (4) 


Second, we notice that in our cubic spline interpolation 
method the slope of the tensor spectrum becomes larger 
in the low-k regime, but still is consistent with the power- 
law parameterization in 20 confidence level. 

As we argued in the previous section, there exists an 
extraordinary power excess in the higher wavenumber 
regimes of BICEP2 data. Given this consideration, in 
what follows we only adopt the selected first 5 band- 
power data of BICEP2 for our reconstruction. Based 
on the multipole ranges covered by these powerbands 
(20 < £ < 200), we sample points at kı = 0.002, 
k2 = 0.0063 and k3 = 0.02 Mpc™!, respectively. 


5 


— ACDM + InB, + InB, + InBy 
=: — ACDM+m+r a 
4b O k=0.002,In[10!°P,(k)] = —0.3602 oo 


3b 


In{10! P;(k)] 


k{Mpe!} 


FIG. 4: Reconstruction of primordial tensor spectrum with 
BICEP2 data (5 bandpowers). 


The reconstructed tensor spectra as well as error bars 
are shown in Fig. First, from the power-law param- 
eterization, we can see that the blue tensor spectra are 
still favored but with the tilt becomes smaller as expected 


ro.002 < 0.067 ,95%CL , (5 bandpowers), (5) 
ne = 1.201986 ,95%CL . (5 bandpowers) (6) 


Second, with only the first 5 bandpower data, unlike the 
simplest power-law parameterization, a non-trivial shape 
of tensor spectrum is obtained. Concretely, in the range 
of k € (0.002,0.006) BICEP2 data favor a large tensor 
blue-tilt, while when k > 0.0063Mpc~! the spectrum be- 
comes almost flat. This is due to the fact that we do 
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not use the last 4 bandpower data. The resulting tensor- 
to-scalar ratio in our cubic spline interpolation method 
is 


0.002 < 0.060 ,95%CL , (5 bandpowers) . (7) 


r r 
—— ACDM + InB, + InBy + InB; : Planck + WP 
4F — — — ACDM + InB, + InB + InB; : Planck +7 prior 


k|Mpe!] 


FIG. 5: Reconstruction of primordial tensor spectrum with 
Planck+WP and Planck+7 prior. 


Because the above blue tensor tilt is significantly incon- 
sistent with the standard inflationary prediction, we have 
to figure out its reason on the data analysis level. No- 
tice that for now we always utilize Planck+WP-+BICEP2 
data compilation, so one natural guess is that this 
blue tensor tilt might reflect the tension among Planck, 
WMAP polarization and BICEP2 data sets. In order to 
justify our conjecture, we have to remove the above data 
sets one by one. Since we vary both the primordial spec- 
trum and standard ACDM parameters, such as baryon 
(Qsh?), cold dark matter density (Nch?) etc, we have to 
keep the robust Planck temperature data in order to get 
well estimations of the ACDM parameters. Hence, we 
first remove the BICEP2, i.e. using Planck+WP, and 
further discard WMAP polarization data, i.e. only using 
Planck temperature data. However, due to the fact that 
CMB temperature data are insensitive to the reionization 
optical depth (7), if we discard WMAP polarization data 
we have to include gaussian prior on 7 to break the well- 
known degeneracy between 7 and the scalar amplitude 
Ag. Here we take the gaussian prior as 


r = 0.089 + 0.013. (8) 


Besides this, because the contribution of tensor spec- 
tra to CMB temperature anisotropies is only significant 
in the mulitpole range (2 < € < 100), when we use 
Planck+WP or Planck+7 prior data sets, we sample the 
k knots of tensor spectra in the range of (0.0002, 0.01). 

The resulting primordial tensor spectrum shape, cor- 
responding marginalized 1D/2D posterior distributions 
and parameter regimes are shown in Fig] Fig[I0] and 
Tab[V] First, we can see that the reconstructed ten- 
sor spectra from Planck+WP and Planck+7 prior are 


r r 
ACDM + InB, + InB, + InB; : Planck + WP + BICEP2(9bins) 
= — —ACDM + InB, + InB, + InBy : Planck + WP 


In{10"" P;(k)] 


k|Mpe!] 


FIG. 6: Reconstruction of primordial tensor spectrum 
with/without BICEP2. 


very similar, both the central values and the marginal- 
ized error bars. So, we can draw the conclusion that 
WMAP low-/ polarization data are not crucial for the 
tensor reconstruction results. Second, as shown in Fig[6] 
the error bars from Planck temperature data are quite 
large compared with those from the data compilation 
Planck+WP+BICEP2 (see the error bars at sampling 
knot k = 0.01Mpc™! in the dashed black curve and those 
at the second knot in the black solid curve in Fig[6). 
It means that the current Planck temperature data are 
not rubost to determine the shape of the primordial ten- 
sor spectrum. The resulting tensor spectrum from only 
Planck temperature data could be any shape among red, 
blue or scale-invariant types. Third, when we compare 
the second knot in Planck+WP (dashed black curve) and 
the first left knot Planck+WP-+BICEP2 results (solid 
black curve), Fig[6] shows that both their central value 
and marginalized error bars are very close. It reflects the 
fact that the reconstructed tensor spectrum in the low- 
k regime is actually driven by Planck temperature data. 
Furthermore, considering the fact that the BICEP2 data 
dominate the high-k part, we conclude that our recon- 
structed blue tensor tilt are due to the tension between 
Planck and BICEP2 data sets, i.e. the fact that small 
tensor amplitude signals from Planck temperature data 
dominates the large scale reconstruction, while the large 
tensor amplitude signals from BICEP2 B-mode polariza- 
tion data dominates the small scale reconstruction, leads 
to the resulting blue tensor tilt. 

Finally, in order to give an intuitive impression of our 
reconstruction result, in Fig. [7] we plot the BB auto- 
correlation power spectrum of our marginalized mean 
models (listed in Tab. of Appendix as well as 
the scale-invariant tensor spectrum with r = 0.2 model 
against the BICEP2 bandpowers data sets. We can see 
that in order to fit the BICEP2 data in the last 4 band- 
powers, compared with the 5 bandpower reconstruction 
result (green curve) and the scale-invariant one (blue), 
the 9 bandpower (red) curve grows up in the high-é 
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cubic spline (BICEP2: 9 bandpowers) 


cubic spline (BICEP2: 5 bandpowers) ------- 
ACDM+r=0.2 se 


BICEP? datasets == 


1+ 1)C8® 2n[ WK? ] 


1 1 1 L 1 1 
50 100 150 200 250 300 350 


FIG. 7: CMB B-mode polarization auto-correlation spectrum 
OF. 


regime significantly. 


IV. CONCLUSIONS 


Starting with a purely phenomenological point of view, 
in this paper we have reconstructed the shape of the 
primordial scalar and tensor spectra by using the cu- 
bic spline interpolation method with Planck tempera- 
ture and BICEP2 B-mode polarization data sets. We 
find that, due to the anti-correlation between scalar 
and tensor amplitudes on the large scales, the large 
value of tensor-to-scalar ratio discovered by BICEP2 data 
will lead to the suppression of scalar amplitude in this 
regime. Concretely, the vanishing scalar index running 
model is strongly disfavored by Planck+WP-+BICEP2 
data compilation with more than 30 confidence level 
on the k = 0.0002 Mpc™! scale. Furthermore, for the 
tensor spectrum reconstruction, a blue-tilt spectrum is 
obtained no matter using only the first 5 bandpow- 
ers ne = 1.207925 (95%CL) or the full 9 bandpowers 
ny, = 1.244928 (95%CL) of BICEP2 data sets. Because 
of the large tensor tilt, compared with the large tensor- 
to-scalar ratio value (r ~ 0.20) under the scale-invariant 
assumption, our cubic spline interpolation method gives 
ro.oo2 < 0.060 (95%CL) and 19.902 < 0.064 (95%CL) by 
using the data sets Planck+WP+BICEP2 (5 bandpow- 
ers) and (9 bandpowers), respectively. Finally, we also 
studied the data without BICEP2, we found that our re- 
sulting blue tensor tilt actually reflects the tension in the 
tensor amplitude between Planck (small amplitude but 
dominate the reconstruction on the large scale) and BI- 
CEP2 (large amplitude but dominate the reconstruction 
on the small scale) data sets. 

Our results show that the conclusion of the blue-tilt 
tensor spectrum is very significant and independent of 
using power-law or cubic spline parameterizations. More 
important, this blue-tilt spectrum is not consistent with 
the prediction of the standard single field inflationary 
paradigm ng = —r/8. On the one hand, it asks for a 


more careful cross-check with future experiments, such 
as the polarization data of Planck and Keck Array. On 
the other hand, once this discovery is confirmed, it will 
lead to a paradigm revolution about our understanding 
of the early universe. 
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Appendix A: Marginalized statistics in scalar 
spectrum reconstruction 


Here we list the various marginalized statistical results 
for cubic spline interpolation and power-law parameter- 
izations of scalar spectrum, including 1D, 2D marginal- 
ized posterior distribution, marginalized mean values as 
well as the 68% (or 95%) confidence levels. 


— ACDM-n,-In Agtr+in A;+In Ag+in Ag: Planck+WP+BICEP2 
— ACDM-ng-In Astr+in Ay+In Ag+in Ag: Planck+WP 
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FIG. 8: 1D/2D posterior distribution of scalar spectrum reconstruction. 
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Parameters mean 4 


2.2244 
0.12044 
1.041254 


< 0.41 ( 


E 68% C.L. 
0.031 
-0.0028 
-0.00065 


95%CL) 


mean a 


t 68% C.L. 


mean + 68% C.L. 


2.2242 
0.12023 
1.041332 


0.030 
-0.0027 
0.00063 
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FIG. 9: 1D/2D posterior distribution of tensor spectrum reconstruction with BICEP2. 
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